ADVANCED VISION INTERVENTION ALOGRITHM 

1 Cross Reference to Related Application : 

2 This application is a continuation of Provisional Patent Application Serial No. 60/250,296, 

3 filed November 30, 2000. 

4 Field of the Invention: 

5 This invention relates in general to the field of optically corrected aberrant human vision and 

6 more particularly, relates to the construct of an improved optical model of the human eye and 

7 procedures to achieve a more precise optical correction. 

ip Background of the Invention : 

pi The human eye often contains aberrations that may be optically corrected for normal 

10 * activities. If the primary aberration causes a distant object to be focused in front of the retina, for 

11* example, the eye is said to be myopic or near-sighted. If the primary aberration causes the distant 

ill object to be focused behind the retina, the eye is said to be hyperopic or far-sighted. Often the 

iT aberrations are such that the light rays which enter the cornea in one location from a distant object 

14 focus at one point inside the eye while light rays which enter the cornea in another location from the 

15 same distant object focus at a different point inside the eye. A typical example of this is called 

1 6 astigmatism in which light rays which enter the eye in one meridian focus at one point and light rays 

1 7 which enter in the meridian perpendicular focus at another location. It is easy to visualize that other 

18 more general aberrations can and do occur in the human eye. 

19 Interventions in common use today to correct such human eye aberrations are spectacles, 
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1 contact tenses, corneal refractive surgery, such as LASIK or corneal implants, and phakic intraocular 

2 lenses. The correct specification of spherocylindrical spectacles and contact lenses is well 

3 established. The newer surgery based techniques are not as predictable as may be desired, but still 

4 are in wide spread use. 

5 A current development is instruments intended to measure the total aberrations of the eye 

6 rather than simply the traditional spectacle correction parameters of sphere, cylinder, and axis. These 

7 instruments provide the theoretical information required to improve on vision correction beyond the 

8 ' lower-order aberrations of sphere, cylinder, and axis. The application of this information in a clinical 

9 setting has not been entirely successful because the optical measurements of the eye and the 
1|| intervention methods are not exact. 

Ml Another area of interest is the relatively mature process of calculating the intraocular lens 

]|* (IOL) power for post cataract patients. The calculation of implantable lens powers and prediction 

l-f of post-operative outcomes has been an ongoing research interest since the first implanted lens. The 

pi calculations employed by modern strategies may effectively reduce the prediction error bias for the 

M majority of normal cases but even with these adaptive third-generation IOL calculations, there are 

ft a significant number of real world cases whose outcomes are poorly predicted. A particular category 

17 that is poorly predicted is cataract surgery following previous refractive surgery. A primary source 

18 of error in the corneal power measurement with keratometers is that these instruments typically 

1 9 measure corneal curvature in the 32 mm zone and can miss the relatively flat central region. Modern 

20 IOL calculation strategy, as stated above, is based on calculations involving a mix of parameters 

2 1 intended to improve the effective lens position predictor (ELP) and calculation. Attempts to improve 

22 ELP have employed such parameters as axial length, eomeal power, horizontal corneal diameter, 
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anterior chamber depth, crystalline lens thickness, preoperative retraction, and age. The basic 
methodology for the modern IOL calculation strategy has three steps: a) compute the ideal theoretical 
power of the implanted lens (this lens power will not necessarily be available) using measured and 
adapted optical model parameters, b) for the surgeon selected available power near the ideal power, 
predict the post operative refractive outcome, and c) given the actual outcome of the procedure, 
update optical model parameters to decrease the prediction error of future surgeries. As stated above, 
such modern strategies, while improving the prediction error bias for the majority of normal cases, 
do not improve predicability for many other cases, such as the case of cataract surgery following 
previous refractive surgery. One modern attempt at improvement is custom LASIK using corneal 
topography and/or wavefront aberration data. However, large amounts of positive spherical 
aberration can be induced as a result of LASIK procedures and for large pupil diameters, the retinal 
image quality may be degraded. 

While past and current developments have improved the predictability of optical correction 
of aberrations in the human eye, as may be understood by those skilled in the art, the presently 
known techniques or methodologies have not produced a measure of predictability for the many 
variations found in the parameters describing the human eye. 

Summary of the Invention : 

The inventive principles of the advanced vision intervention algorithm (A VIA), as disclosed 
in the description of the preferred embodiment, may be used to improve the predictability of most 
known methodolgies for optical vision correction and may be applied to any similar methodology 
which may be developed, The advanced vision intervention algorithm (A VIA) provides a means to 



optimize the predictability of almost any current or anticipated customized human vision 
intervention method. In the following summary and Description of a Preferred Embodiment, AVIA 
is described in general terms in a series of steps. For the purpose of explanation, the AVIA method 
is shown in nine steps but as would be known by those skilled in the art, the number of steps used 
to describe the invention may be varied without departing from the principles of the disclosed 
invention. 

For the example shown in the preferred embodiment, a summary of these steps is: 

1. Input data: The input data is evaluated for obvious errors. Individual classes are used 
for each data type so that new data types and their evaluation methods are easily added to the 
framework. This includes inputting a set of categorical data for an eye and visual correction method 
and inputting a set of continuous data for an eye and visual correction method. The input data 
expresses the patient's visual correction need in terms of wavefront aberration. If the wavefront 
aberration is not explicitly measured using a system intended for this purpose, it may be estimated 
using the information provided such as spectacle 

2. Build Optical Model: The second step is building an optical model of the eye and 
visual correction method that has been adapted based on prior visual correction outcomes. This 
processing step employs a ray transfer element (RTE) that is described in detail below. The RTE is 
an important component of AVIA. 

3. Generate Initial intervention recommendation. In this step an initial intervention 
recommendation is computed based upon the optical model. This recommendation will either be 
categorical (e.g., selection from a series of available implantable lenses); a continuous scalar such 
as a theoretical power value; or a continuous vector quantity such as a surface description, or a 



1 combination of these such as a lens series with a discrete set of base curves and a customizable front 

2 surface. 

3 4. Outcome Prediction: In this step, the outcome of the intervention is predicted in terms 

4 of the wavefront aberration, The predicted wavefront aberration is displayed in a number of ways 

5 any one of which can be selected by the surgeon for viewing. Examples of these displays include 

6 equivalent spectacle correction, wavefront aberration variance, wavefront aberration contour map, 

7 point spread function metrics and displays, modulation transfer function metrics and displays, and 

8 simulated retinal images. 

9 5. Predicted Outcome Evaluation: In this step the surgeon uses the various displays of 
lil the predict ed outcome wavefront aberration to decide if the outcome would be acceptable for the 
ljf- current case under consideration. 

6. Intervention Design Iteration: In this step, the surgeon may continue iterating a 

13 modification of the Intervention plan and re-evaluates the predicted outcome until the predicted 

lf# outcome is judged as acceptable. 

15: 7. Perform Intervention: In this step the surgeon performs the intervention. 

fijf 8. Evaluation of Outcome: In this step the outcome of the intervention is assessed using 

17 data such as wavefront aberration exams, uncorrected and best corrected visual acuity, contrast 

18 sensitivity test, post-intervention spectaeie correction, etc. 

19 9. Update Historical Database: In this step the database of historical visual correction 

20 cases is updated to include the current case. 
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1 Brief Description Of The Drawings : 

2 Figure 1 shows in a block element flow diagram, the method of a preferred embodiment, 

3 according to the disclosed inventive principles; 

4 Figure 2 shows in a simplified optical system for the human eye, how step 2, as described 

5 in the Summary and Description of a Preferred Embodiment and shown as numeral 2, in Figure 1 , 

6 for the modeling and use of a RTE; 

7 Figure 3 shows an intersection point of an incident ray 1, and the RTE referenced to the 

8 meridional and transverse planes; and 

9 Figure 4 shows the basis vectors for the meridional coordinate system of Figure 3. 

\4i 

li s ; Description of a Preferred Embodiment : 

The invention, according to its inventive principles, is disclosed with reference to the 

1=3* preferred embodiment, as described below, and with reference to the accompanying drawings. As 

ll# shown in step 1 in Figure 1 , the first step uses the clinical and categorical data of the eye to construct 

15 : an optical model, as shown in Figure 2, and in step 2 of Figure 1 . 

16 Optical Model 

17 The foundation of the optical calculations and predicted outcome of a given intervention 

1 8 strategy is the optical model. The optical model is created using all available data for a given eye. 

19 This data could be as simple as the two basic quantities of keratometric data and spectacle correction 

20 or as complex as specifying multiple surfaces, axial locations, and wavefront aberration data, or 

21 additional input data, without departing from the disclosed inventive principles. To explain the 

22 procedure, a simple model consisting of a corneal surface and wavefront aberrations only, is used. 
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1 In Figure 2 is shown the cornea and wavefront W positioned along an optical axis. In this Figure 2, 

2 the wavefront is located at the entrance pupil of the eye. As shown in Figure 2, if a ray R2 is traced 

3 from the fovea, it would appear to exit the eye in the direction of RO. This exiting ray is 

4 perpendicular to W as shown by the dotted line intersection with W. 

5 To construct a model consistent with this information an average axial length of 24 mm may 

6 be assumed and a ray transfer element (RTE) introduced that maps an incoming ray R2 from the 

7 fovea to the intermediate ray Rl so that the exit ray RO has the desired orthogonality relationship 

8 with W, as shown in Figure 2. 

9 Construction of the RTE is accomplished using the following steps: 
10; 1 . For a regular sampling of points on W do the following steps: 

1 C 2. Construct the ray RO perpendicular to W at the sample point on W; 

ftl 3 . Trace - RO through the system until it intersects the RTE plane; 

if* 1 4. At the RTE intersection point, compute the mapping for the incident and transmitted rays; 

|fc and 

|g 5. If there are more samples on W go back to step 2. 

f# The optical ray tracing of refractive and reflective elements is well known to those skilled 

17 in the art, for example as shown in Modern Optical Engineering by W.J. Smith and published by 

1 8 McGraw-Hill, 1 990, and for that reason is not described in detail. 

19 Ray Transfer Element 

20 The ray transfer element (RTE) is an important aspect of AVIA. It provides both a means to 

21 generate an optical model that is consistent with exam data and a means to adapt the modeling to 

22 past surgical procedures so that prediction of postoperative outcomes is improved for future 
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surgeries. The basic operation of the RTE is to efficiently map incident rays to transferred rays 
(either transmitted or reflected rays) and provide a method to smoothly adjust the mapping between 
data points to simulate continuous data. Primary operations on the RTE include ray intersection and 
ray reflection and transmission. The physical structure of the RTE is modeled as a plane with an 
aperture. Usually the aperture is circular, so the RTE looks like a circle of zero thickness in space. 
In local coordinates the RTE is represented as a plane equation with a separate test for the aperture. 
The equation of the planar object In local coordinates is simply 

Z = 0 

(1) 

The intersection of the RTE with a ray is straightforward. Let a ray in local coordinates be defined 
as In (2). 



X 




Xo 




Xd 


y 




yo 


+ t 


yd 


z 




Zo 




Zd 



Then the interseciion of the ray with the plane may be found, as in equation (3). 
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where z d = 0. Given the location of the planar Intersection (x,y, 0), a final test is made to ensure the 
point is inside the aperture region. 
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1 RTE Construction 

2 In the preferred embodiment, the construction of the RTE is described in local coordinates. In 

3 practice, rays defined in world coordinates are first transformed into local coordinates using 

4 homogeneous transformation matrices, collectively called a world-to local transformation matrix. 

5 Likewise, resulting rays in local coordinates are transformed to world coordinates using an inverse 

6 of the previous homogeneous transformation matrix, called local-to-world transformation matrix, 

7 as would be well known to those skilled in the art (See for example, Glasser A.S., An Introduction 

8 to Ray Tracing, Academic Press, 1 989, for a discussion of homogeneous transformation matrices.) 

9: Since the goal of the RTE is to transform an input ray into an output ray for a given intersection 

1 0" point, we simply construct a mapping for a set of input/output ray pairs. Suppose an incident ray I 

ijjj and a transmitted ray r defined in local Coordinates is given. The intersection point (x,y, 0) is denoted 

K S. The plane that contains the optical axis (z-axis) and the point S is called the meridional plane. The 

it plane parallel to the z-axis and perpendicular to the meridion at plane is referred to as the transverse 

1 4 : plane. These planes are illustrated in Figure 3 . The strategy for finding the mapping of I to T is to 

i 5 find the change in elevation and azimuth angles for the two rays with respect to the meridional and 

16 transverse pianos, so that I can be rotated into T. The first step is to compute a meridional plan 

17 coordinate basis matrix. For this coordinate system, X points from S to the origin. Y is 90 deg 

1 8 counter clockwise to X, and T is parallel to the original z-axis. These basis vectors are illustrated 

19 in Figure 4. 

20 For simplicity of presentation, in the following the prime notation is dropped for the 

21 meridional basis vectors X', Y, and Z'. Given S we can compute the meridional coordinate basis 
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1 matrix by finding the X, Y, and Z vectors as in equation (4). 
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If the norm of S is zero, the intersection point is at the origin, i.e., on the optical axis. In this case X 
in equation (4) is taken as being along the x-axis. Using this the meridional coordinate basis matrix 
M can be written as In (5). 
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(5) 



To compute the desired rotation angles between / and T, first rotate the vectors using the meridional 
coordinate basis Mas in (6). 



ft 

14 
15 
16 
17 
18 
19 



T M = M r T 



(6) 



In (6) and in the following discussion, the vectors / and T are first normalized to unit length. Now, 
can be computed the rotation angle for each of these vectors to the meridional coordinate x-axis. 
Also can be computed the rotation angle between the XZ plane and each of these vectors. These 
calculations are made using equations (7) and (8), respectively. 
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0a = tan-* 
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In equations (7) and (8), A stands for either the vector / or T. The arc tangent function 
indicated in these equations is conveniently computed using the a atan2 function in C++ to 
automatically handle all possible values of numerator and denominator around the unit circle 
including the case of denominator equal to zero. After these two elevation angles and two azimuth 
angles have been computed, the difference can be computer as indicated in (9). These are the angles 
stored to map an incident ray direction to the transferred ray direction. 



d9 = 6 j — 9 j 

(9) 

da - a j- a T 

To summarize, for a given intersection point S, compute and save the elevation and azimuth 
angles d6 and da and do required to rotate the incident vector I into the direction of transferred vector 
T. All of the forgoing calculations would be known to those skilled in the art and for that reason are 
not described in detail. Discussed next, is the use of the RTE to map an input vector to an output 
vector. This illustrates the case of forward ray tracing. The case of reverse ray tracing is similarly 
computed, as would be known to one skilled in the art and is omitted for brevity. 
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1 RTE Optical Ray Tracing 

2 To perform an optical ray tracing of the RTE, first find where an incident vector I is 

3 intersected with the RTE at location S. This intersection point is computed using equation (3). Given 

4 this intersection point, the transmitted vector T is easily computed. In general, the intersection point 

5 S will not correspond to a previous constraint point, so some type of interpolation will be required. 

6 Two strategies to accomplish this are (1 ) to compute a least squares approximation function fit (e.g., 

7 2D polynomial or 2D B-Spline) for each of the azimuth and elevation angles or (2) to perform a 

8 Delaunay triangulation of the intersection points in the XY plans and interpolate the points to obtain 

9 the angles. As would be known to those skilled in the art.., Suitable software may be developed for 
1 0= performing these calculations on a general purpose computer. The computation of 7/(which produces 
] f a ray tracing of the RTE) is performed as follows: 

E2 1. Intersect I with the RTE to find S 

{f: 2. Apply the transformations indicated in (4) to (6) and then rotate I by the azimuth and 

% elevation angles associated with S. This rotated vector represents T in meridional coordinates. 

15 3. Apply the inverse coordinate transformation to T. 

16 The adjustment angles dB and da are applied to the angles interpolated in step 2 of the 

17 algorithm. This provides a method to adjust the optical model based on actual outcomes. In a similar 

18 way, a paraxial power value could be assigned to the projection of the I end T rays onto the 

19 meridional and transverse planes. This alternate scheme has the benefit of describing the bending 

20 of rays in units of diopters that are familiar to visual optics professionals. 
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1 Strategy for RTE Updating Based on Historical Cases 

2 A strategy for updating the RTE portion of the optical model is based on the evaluation of 

3 historical cases. In a preferred embodiment, as described in this application, first examine the 

4 historical database for similar cases to a given present case under consideration. This similarity is 

5 measured in terms of both categorical data and continuous data. Categorical data such as surgeon, 

6 sex, intervention procedure, etc., should be matched closely by all historical cases to be similar to 

7 the current case. This is accomplished by filtering the database for all records that match the same 

8 categories. Continuous data such as preoperative spectacle correction, corneal thickness, etc., is 

9 measured using a weighted Euclidean distance. Each of the 5 continuous data is normalized so that 
Y§ the mean and standard deviation of each random variable is zero and one, respectively . Only the N 
1 f closest cases from the historical database are used to compute the RTE adaptation parameters. A 
]|l reasonable value for the integer N may be determined during a first evaluation (described below). 

it The RTE update algorithm is then: 

M 1 . Establish the categorical data for the current case under consideration and identify 

15 the continuous data. 

16 2. Filter the historical database so that we only consider cases in the same category. 

17 3. For all filtered records, find the N closest records using normalized distances. 

18 4. For these N closest records, find the average back-calculated RTE for the actual 

19 postoperative wavefront aberration (See RTE Arithmetic Operations below.) 

20 5 . Use the average of the RTE deltas in the Optical Modeling of the current case 
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1 I It is often possible to have the case of time-varying parameters, that is, cases where the 

2 statistics slowly vary over time. To account for this in our calculations of the updated RTE, a 

3 weighting may be applied to the average calculation ab that older cases are not weighted as much 

4 as more recent cases. 



5 RTE Arithmetic Operations 

6 RTE arithmetic operations that help facilitate AVIA calculations are: 

7 Addition, subtraction, multiplication, and division of two RTEs. 

8 Addition, subtraction, multiplication, and division of an RTE and a scalar. 

|f These operations (addition, subtraction, multiplication, and division) on two RTEs are 

ill performed in a point-wise fashion. In general, the stored sample points from two RTEs will not be 

10] in the same locations, so one or both of the sample point arrays must be interpolated to obtain points 

n at matching locations. 

if- These operations (addition, subtraction, multiplication, and division) on a RTE and a scalar 

H are also performed in a point-wise fashion. In this case there is no need to interpolate to align sample 

15 points as the same scalar operation is applied to all points in the RTE. 

16 Implementation of the RTE using the C++ programming language permits these operations 

17 to be specified directly using standard symbols of "+", "- . and "F. This makes the generation and 

18 maintenance of RTE aware programs very efficient. For example, the mean of four RTEs can be 

19 expressed as: 
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(R1+R2+R3+R4) 

MR = (10) 
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Using an appropriate class implementation of a RTE, this expression is a valid C++ statement 
and its intent is easily understood. 
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